Topology hidden behind the breakdown of adiabaticity 
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For classical Hamiltonian systems, the adiabatic condition may fail at some critical points. How- 
ever, the breakdown of the adiabatic condition does not always make the adiabatic evolution be 
destroyed. In this paper, we suggest a supplemental condition of the adiabatic evolution for the 
fixed points of classical Hamiltonian systems when the adiabatic condition breaks down at the crit- 
ical points. As an example, we investigate the adiabatic evolution of the fixed points of a classical 
Hamiltonian system which has a number of applications. 

PACS numbers: 05.45.-a, 03.75.Fi, 02.40.Pc 



I. INTRODUCTION 

Adiabaticity is an interesting concept in physics both 
for theoretical studies and experimental practices 0, 0, 
0,0,0. According to the adiabatic theorem Q], if the pa- 
rameters of the system vary with time much more slowly 
than the intrinsic motion of the system, the system will 
undergo the adiabatic evolution. For a classical system, 
the adiabatic evolution means that the action of the tra- 
jectory keeps invariant. For a quantum system, an ini- 
tial nondcgenerate eigenstate remains to be an instanta- 
neous eigenstate when the Hamiltonian changes slowly 
compared to the level spacings Hence, the adiabatic 
evolution has been employed as an important method of 
preparation and control of quantum states 0, Q, S 13 • 

However, a problem may arise when the eigenstates 
become accident degenerate at a critical point, i.e., when 
the level spacing tends to zero at a critical point. For a 
classical system it corresponds to that the frequency of 
the fixed point is zero at the critical point. The adiabatic 
condition is not satisfied at the critical point because the 
typical time of the intrinsic motion of the system becomes 
infinite. Can adiabatic evolution still hold up when the 
adiabatic condition breaks down at the critical point? 

Our motivation, derives from practical applications in 
current pursuits of adiabatic control of Bose Einstein 
condensates (BECs) ^j, which can often be accurately 
described by the nonlinear Schrodinger equation. Here 
the nonlinearity is from a mean field treatment of the 
interactions between atoms. Difficulties arise not only 
from the lack of unitarity in the evolution of the states 
but also from the absence of the superposition princi- 
ple. This was recently addressed for BECs in some spe- 
cific cases But then, however, for such systems, 
only finite number of levels are concerned. The nonlinear 
Schrodinger equation of the system with finite number of 
levels can be translated into a mathematically equivalent 
classical Hamiltonian system. The evolution of an eigen- 
state just corresponds to the evolution of a fixed point 
of the classical Hamiltonian system. Then, the accident 
degeneracy of eigenstates is just translated into accident 
collision of the fixed points. The latter one is quite well- 



known subject and has been studied widely at least as a 
purely mathematical problem |l3j|. Hence, our concern 
here is only focused on the adiabatic evolution of the 
fixed points of classical Hamiltonian systems. 

In this paper, we present a supplemental condition 
of the adiabatic evolution for the fixed points of clas- 
sical Hamiltonian systems when the adiabatic condition 
breaks down at some critical points in the terms of topol- 
ogy. As an example, we investigate the adiabatic evolu- 
tion of the fixed points of a classical Hamiltonian system 
which has a number of practical interests. We show that 
the adiabatic condition will break down at bifurcation 
points of the fixed points. But the adiabatic evolution is 
destroyed only for the limit point. For the branch pro- 
cess, the adiabatic evolution will hold, and the correc- 
tions to the adiabatic approximation tend to zero with a 
power law of the sweeping rate. 



II. SUPPLEMENTAL ADIABATIC CONDITION 
FOR THE FIXED POINTS OF CLASSICAL 
HAMILTONIAN SYSTEMS 

For clarity and simplicity, we consider a one-freedom 
classical Hamiltonian H (p, q; A) with canonically conju- 
gate coordinates (p, q) where A is a parameter of this 
system. The equations of motion are: 



dH 

dp ' 



dH 



(1) 



We can find two kinds of trajectories in the phase space 
for the system: fixed points and closed orbits. The fixed 
points are the solutions of Eqs. Q when the right hands 
of them are zero. For a Hamiltonian system there are 
only two kinds of the fixed points: elliptic points (stable 
fixed points), hyperbolic points (unstable fixed points). 
The closed orbits are around each of the elliptic points. 
We denote the fixed points by z*(p,q) (i — 1,2, ••• , I) 
where I is the total number of the fixed points. 
The action of a trajectory is defined as 



1 = 2tt ^ ! " ,q ' 



(2) 
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where the integral is along the closed orbit. Obviously, 
the action of a fixed point is zero. The action is invariant 
when system undergos adiabatic evolution. 

According to the adiabatic theorem pj, the adiabatic 
condition can be expressed as 



2tt dX 

Ultt << ' 



(3) 



where is the frequency of the fixed point. If this condi- 
tion holds, the system will undergo adiabatic evolution, 
and keep the action not varying. If O ^ 0, the condition 
can always be satisfied. 

We can obtain the frequencies of the fixed points by 
linearized the equations of motion. Let us define the 
Jacobian matrix as 



J 



d 2 H d 2 H 

da' 2 dqdp 

d 2 H d^H 

dqdp dp 2 



(4) 



It is well-known that when det(J)| z * > 0, the fixed point 
is a stable fixed point (elliptic point); when det( J)| z * < 0, 
the fixed point z* is a unstable fixed point (hyperbolic 
point). The point with det( J)\ z * = is a degenerate point 
at which the stability of the system is not determined. 

For a stable fixed point z*, the frequency of this fixed 
point is 



Oo = V|det(J)| a 



(5) 



Obviously, Slo depends on the parameter A. 

Supposing at a critical point, namely A = A c , we have 
^o(A c ) = 0. Therefore, the condition will break down 
at the point. We want to know what will happen when 
the adiabatic condition fails (will the adiabatic evolution 
of the fixed point be destroyed when the adiabatic con- 
dition does not hold ?). 

In fact, if |det(J)| z « = 0, the point 2* is a bifurca- 
tion point at which the fixed point will collide with the 
other fixed points [ijj [l4(. Hence, the breakdown of 
adiabatic condition is equivalent to collision of the fixed 
points (equivalent to accident degeneracy of eigenstates 
of a corresponding quantum system). In the collision 
process, fixed points may annihilate or merge into a sta- 
ble fixed point. The collision of the fixed points can be 
described clearly in the terminology of topology . 

The equations of motion just define a tangent vector 
field (f>(p, q) — (^jp-, — ^) on the phase space. Obviously, 
the fixed points z*(i = 1, 2, • ■ • I) are the zero points of 
the vector field, i.e., <fr(z*) = 0. We know that the sum of 
the topological indices of the zero points of the tangent 
vector field is the Euler number of the phase space which 
is a topological invariant jTsl ] . For a Hamiltonian system, 
the topological index for a stable fixed point is +1 and 
for a unstable fixed point is —1. 

Indeed, if the fixed point is a regular point (not a de- 
generated point), i.e., det(J)| z * 7^ 0, the topological in- 
dex of the fixed point can be determined by determinant 
of the Jacobian matrix defined by Eq. Q _2 El • If 



det(J)\ z * > 0, z* is a stable fixed point and the topolog- 
ical index is +1; if det(J)| z » < 0, it is a unstable fixed 
point and the index is —1. 

If det(J)| z * = 0, i.e. if z* is a bifurcation point, the 
topological index of this point seems to be not deter- 
mined. As we have shown before, the point is just the 
critical point of adiabatic evolution, corresponding to col- 
lision of the fixed points. 

However, because the sum of the topological indices 
is a topological invariant, the topological index is con- 
served in a collision process of the fixed points. There- 
fore, the topological index of the bifurcation point can be 
determined by the sum of the indices of the fixed points 
involved in collision. So, if the topological index of the 
bifurcation point is not zero, it is still a fixed point after 
collision. But if the topological index of the bifurcation 
point is zero, the bifurcation point will not be a fixed 
point after collision. 

Now, let us imagine what will happen when a fixed 
point is destroyed by a collision process. Because there 
are only two kinds of trajectories for a classical Hamilto- 
nian system: fixed points and closed orbits around each of 
the stable fixed point, so when a fixed point is destroyed, 
it will form a closed orbit around the nearest stable fixed 
point. The action of the new orbit must be proportional 
to the distance between the critical point and the nearest 
stable fixed point. This sudden change of action (from 
zero to finite) is so-called "adiabatic tunneling probabil- 
ity" which has been studied in Refs. |16l Il7j. On the 
other hand, if the topological index of the bifurcation 
point is —1, i.e., it is a unstable fixed point after the 
collision, we can not expect the adiabatic evolution can 
keep on after collision. 

But if the topological index of the bifurcation point is 
+1, i.e., it is still a stable fixed point after the collision, 
or in other word, the stable fixed point survive after col- 
lision. For such case, the adiabatic evolution will not be 
destroyed. 

From above discussion, it is clear that when the adi- 
abatic condition given by Eq. Q does not hold at a 
critical point with Q,q = ^/| det(J)| z *(^ c ) = 0, the system 
will still undergo the adiabatic evolution if the topologi- 
cal index of the fixed point z*(A c ) is +1. On the contrary, 
if the topological index of the point z*(A c ) is zero or —1 
the adiabatic evolution will be destroyed. 

Hence, we get a supplemental condition of the adia- 
batic evolution of the fixed points for a classical Hamil- 
tonian system when the adiabatic condition breaks down 
at a critical point. When the adiabatic condition is not 
satisfied at a critical point, the topological property of 
the bifurcation point plays an important role to judge 
whether the system will undergo adiabatic evolution over 
this critical point: if the index of the degenerated point 
z*(A c ) is +1 the adiabatic evolution will hold. If the index 
of the point z»(A c ) is zero or —1, the adiabatic evolution 
will not hold. 
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III. A PARADIGMATIC EXAMPLE AND 
APPLICATION 

As a paradigmatic example, we consider the following 
system 

H(z,e,\,-f) = -VI - z 2 cos(9- ^z 2 + jz, (6) 

in which (z, 9) are canonically conjugate coordinates, and 
A, 7 are two parameters. The equations of motion are : 
0=^-, z = --If, these yield 

9 = ^-L_ cos8- Az + 7, (7) 



i = -Vl - z 2 sin 9. (8) 

This classical system can be obtained from a quantum 
nonlinear two-level system, which may arise in a mean- 
field treatment of a many-body system where the parti- 
cles predominantly occupy two energy levels. For exam- 
ple, this model arises in the study of the motion of a small 
polaron Ha. a Bose-Einstein condensate in a double-well 
potential |l9l l2fj l2ll | or in an optical latti ce jg jj 12^ , or 
for two coupled Bose-Einstein condensates j2J, |25j , or for 
a small capacitance Josephon junction where the charg- 
ing energy may be important. This quantum nonlinear 
two-level model has also been used to investigate the spin 
tunneling phenomena recently pr| . 

The fixed points of the classical Hamiltonian system 
are given by the following equations 

6>*=0,7T, 7 - \Z* + — Z rnafl' = (9) 

yl — z* 2 

The number of the fixed points depends on the nonlin- 
ear parameter A. For weak nonlinearity, A < 1, there 
exist only two fixed points, corresponding to the local 
extreme points of the classical Hamiltonian. They are 
elliptic points located on lines 9* — and ir respectively, 
each being surrounded by closed orbits. For strong non- 
linearity, A > 1, two more fixed points appear on the 
line 9* = in the windows — 7 C < 7 < j c , one is ellip- 
tic and the other one is hyperbolic as a saddle point of 
the classical Hamiltonian, where j c = (A 2 / 3 — l) 3 / 2 . In 
the following, we only consider the cases in the region 
-7c < 7 < 7. 

We can obtain the frequencies of the fixed points by 
linearized the Eqs. JJJ and JSJl. For the elliptic fixed 
points on line 9* = 0, the frequencies are equal, they are 

n = y/i - A(i - z* 2 ) 3 / 2 . (10) 



Obviously, if z* — z* — yl — (j) , the frequencies 
will be zero, i.e., flo(z*) — 0. From Eq. 10, we can 



obtain, when A = A c = (7 2// ' 3 + 1) 3//2 , one of the elliptic 
fixed point will be (z* = z*,9* = 0). At this point the 
adiabatic condition will break down. 

Hence, if we start from this elliptic fixed point on line 
9* = at A = Ao > 1, and A changes with time as A = 
Xo — at (keeping 7 invariant in the window — 7 c (Ao) < 7 < 
7c(Ao)), the adiabatic condition © will break down at 
the point z* when A reaches A c because flo = 0. We want 
to know what will happen when the adiabatic condition 
is not satisfied (will the adiabatic evolution be destroyed 
when the adiabatic condition does not hold ?). There are 
two different cases for discussing: 7^0 and 7 = 0. 
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FIG. 1: Time evolution of z(t) initially on the elliptic fixed 
point z*(Xo). (a) for 7 = 0.2, Ao = 2.0. The solid line denotes 
the time evolution of z(t) for a = 0.001 and dotted line for 
a = 0.00001. The dashed line refers to the trajectory of 
the hyperbolic point, (b) for 7 = 0, Ao = 2.0. The solid 
line denotes the time evolution of z(t) for a — 0.000001. The 
dotted line refers to the trajectory of the hyperbolic point. 
The dashed line refers to the trajectory of another elliptic 
point. 

Case 1 (7 ^ 0): For the convenience, we choose 
7 = 0.2 and Ao = 2. We start at the elliptic fixed point 
(z* = z*(Xq),9* — 0) and A varies with very small a. 
At the beginning, the system follows the z*[A(i)] adia- 
batically. But when A reaches A c = (7 2 ^ 3 + l) 3 ' 2 , the 
adiabatic evolution is destroyed with a jump of action 
(the action changes to a finite value from zero suddenly) 
at the point z* (w 0.5048). Fig. 1(a) shows this process. 
Obviously, the breakdown of adiabatical condition leads 
to the destroy of the adiabatic evolution. 

Case 2 (7 = 0): From Eq. |J5J and @, we can have 
two elliptic fixed points on line 9* = for A > 1, 

4 = ±v/l-l/A 2 , 9* = 0, fi = \J A(A 2 — 1), (11) 

and for A < 1, there is only one fixed point, 

z* = 0, 9* = 0, Q = VI - A. (12) 
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Obviously, for z* = and A c = 1., fio = 0, so adiabatic 
condition can not be satisfied. We integrate the classical 
equations of the Hamiltonian system JHJ , with the initial 
condition A = 2, z = z* + (A ), and 9*(0) = 0. Fig. 1(b) 
shows the time evolution of this fixed point for a very 
small sweeping rate a. The final state is a very small 
oscillation around the fixed point (z* = 0,8* = 0). In 
Fig. 2, we plot the dependency of the small oscillation 
amplitude S on the sweeping rate a. From this figure, it 
is clear to see that the amplitude of the small oscillation 
will tend to zero with the sweeping rate decreasing as a 
power law: S = 0.7 * a 1 / 2 . Therefore, for this case, the 
system will evolve adiabatically and keep the action not 
changing for all the time if the sweeping rate is small 
enough, even when A crosses the critical point A c = 1 , at 
which Qq(z*) = 0, i.e., though the adiabatic condition is 
not satisfied when A crosses the point A c = 1, the system 
is still undergoing adiabatic evolution. 



Then, let /j, vary very slowly. At the critical point /i = 0, 
the two wells merge into a single well. At this time, the 
bifurcation point is (x,p) = (0, 0), which is the bottom of 
the single well. So if /i varies very slowly, one can imagine 
that the particle will stay at the bottom of the well all 
the time, even when the system goes from a double well 
to a single one (at this time the adiabatic condition does 
not hold but the bifurcation point is still a stable fixed 
point, because the bifurcation point still corresponds to 
the bottom of the well) . 

As we have discussed in Sect. II, the breakdown of the 
adiabatic condition (f^o = 0) corresponds to the trajec- 
tory bifurcation, i.e., the points z* is just a bifurcation 
point of the fixed points. The properties of the fixed 
points are determined by the following Jacobian 




(13) 



(z*,0) 




Where 



(z,U) = (^f,-^f)- If the Jacobian D ^ 0, 
the zero point (fixed point) is a regular point. But when 
D = 0, the zero point is a bifurcation point. 

There are two kinds of bifurcation points: limit points 
and branch points. The limit point satisfies that D = 

but /J 1 = det 0£ 

V do de / ( a . i0 ) 
sponds to generation and annihilation of the fixed points. 



FIG. 2: The dependency of the oscillation amplitude of the 
final state on the sweeping rate a for the Case 2. The initial 
state and parameters are the same as the Fig 1 (b) except 
the sweeping rate. The solid quadrangles are our numerical 
calculations. The solid line of insert figure is the function 
0.7*a 1/2 . 

In fact, if we make the series expansion of the Hamil- 
tonian © around the critical point, th e sy stem can be 
approximated to a double well system [2(J]. Therefore, 
the phenomenon of Case 2 can be illustrated by the 
standard double well model. Considering a particle in a 
double well, the system is described by the Hamiltonian 
H = 1 /2p 2 - 1 /2^x 2 + 1 /4a; 4 . For fi > 0, it has two stable 
fixed points (x,p) = (•>//!, 0) and (x,p) = (— s/Jl, 0), and 
an unstable fixed point (x.p) = (0,0); for fj, < it has 
a single stable fixed point (x,p) = (0,0). At the criti- 
cal point fi — 0, three fixed points merge into a stable 
fixed point. As the parameter /j, varies from +1 to —1 
the system goes from a double well to a single well. The 
stable fixed points are just the bottom of the wells, and 
the unstable point is just the saddle point of double well. 
If the particle is at the fixed point (y 7 /^ 0) a t the begin- 
ning, i.e., the particle stays at the bottom of one well. 



=/= 0, which corre- 



If D\ 



(z* c ,0) 



D 1 



= 0, the point (z*, 0) is a branch 



point. The branch point corresponds to branch process 
of the fixed points. The directions of all branch curves 
are determined by the equations |14| 



, d 2 z dz „ 

A de +2B dc + c = °> 



c 



dz 2 



2B^+A = 0, 

dz 



(14) 



(15) 



where A, B, and C are three constants. C corresponds 
to A or 7 respectively. Different solutions of the above 
equations correspond to different branch processes. 

For the zero point (z*,0), i.e., the fixed point on line 
6* = 0, we can obtain D = -yi - z* 2 il§| (2 , Q . . Obvi- 
ously, when flo = 0, D = 0, the critical point (z*, 0) is a 
bifurcation point, at which the adiabatic condition fails. 
For the case 1: We can find at the point z*, the 

D 1 = 



Jacobian D = —yl — z 2 V& 



= 0, but 



— syl — z 2 1 / , ^ ^ 0. This point is a limit point which 
corresponds to annihilating of zero points. At this point, 
the elliptic point annihilates simultaneously with a hy- 
perbolic point. In Fig. 1(a), the dashed lines is the tra- 
jectory of the hyperbolic point. Apparently, the elliptic 
point evolves adiabatically until it annihilates with the 
hyperbolic point at z*. After this annihilation the ellip- 
tic point turns to an ordinary closed orbit with a nonzero 
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action, so the adiabatic evolution is destroyed. The anni- 
hilation process of the fixed points of the system JBJ has 
also been discussed in Ref. [22} hi detail. 

For the case 2: At the point z* = 0, the Jacobian de- 
terminant D = 0, and D 1 = —z\J\ — z 2 \, . = 0. This 
is a branch process of the fixed points. We can prove that 
for this case A = C = 0, so the solutions of equations 114|) 
and (JT5j) give two directions: |f = 0, and || = 0. The 
branch process corresponds to merging process. At this 
branch point, three fixed points, two elliptic points and 
one hyperbolic point, merge together. One can see this 
point in Fig. 1(b), in which the dotted line is the trajec- 
tory of the hyperbolic point, and the dashed line corre- 
sponds to the trajectory of another elliptic point. Since 
the total topological index is invariant, the three fixed 
points merge to one point with index +1, i.e., merge to 
an elliptic point. The elliptic point evolves adiabatically 
until it reaches the critical point z* , at which three fixed 
points merges to one elliptic point. Therefore, after the 
branch process, the elliptic point turns to a new elliptic 
point, the action keeps zero and the adiabatic evolution 
still holds. 

From above discussion, we see that the adiabaticity 
breaks down at bifurcation points of the fixed points, 
but only for the limit point the adiabatic evolution is 
destroyed (case 1), while for this case the two fixed points 
annihilate. For case 2, three fixed points merge to one, 
because the critical point (z*,0) is still a stable fixed 
point, the adiabatic evolution keeps with action zero. 

The phenomena discussed above can occur for the adi- 
abatic change of 7 with A fixed. On the other hand, the 
Hamiltonian © is invariant under the transformations 
A — > —A, 9 — > 9 + 7r, and t — > —t. Hence, the phenomena 
can also be found under such transformations. 



IV. CONCLUSION 

In summary, at some critical points, the adiabatic con- 
dition fails, but the adiabatic evolution may not always 
be broken. We find that the topological property of the 
critical point plays an important role for adiabatic evo- 
lution of the fixed points when the adiabatic condition 



does not hold. If the topological index of the critical 
point is +1 the adiabatic evolution of fixed point will not 
be destroyed. On the contrary, if the index of the crit- 
ical point is zero or —1, the adiabatic evolution will be 
destroyed. As a paradigmatic example, we investigated 
the adiabatic evolution of a classical Hamiltonian system 
which has a number of practical interests. For this sys- 
tem, the adiabaticity breaks down at bifurcation points 
of the fixed points. But only for the limit point the adia- 
batic evolution is destroyed. For the branch process, the 
adiabatic evolution will hold, and the corrections to the 
adiabatic approximation tend to zero with a power law 
of the sweeping rate. 

In general, the corrections to the adiabatic approxi- 
mation are exponentially small in the adiabaticity pa- 
rameter, both for quantum system and classical system 
13 It is particularly interesting that the correc- 
tions of the adiabatic approximation may be a power law 
(e.g., for the case 2). The power law corrections to daid- 
abatic approximation have also been found in the non- 
linear Landau-Zener tunneling |16| . In Ref. |l(| . the au- 
thors found that when the nonlinear parameter is smaller 
than a critical value, the adiabatic corrections are expo- 
nentially small in the adiabatic parameter, but when the 
nonlinear parameter equals to the critical value, the adia- 
batic corrections are a power law of the adiabatic param- 
eter. Furthermore, if the nonlinear parameter is larger 
than the critical value, the so-called non-zero adiabatic 
tunneling will occur [161,113. Indeed, the cases, for which 
the corrects to the adiabatic approximation are not ex- 
ponential law with the adiabatic parameter, correspond 
to the collision of fixed points. 
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